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The usual Einstein's equations is modified as a one parameter family of equations in the framework 
of rainbow gravity. In this paper we derive the modified Friedmann-Robertson- Walker (FRW) 
equations when the cosmological evolution of radiation particles is taken into account. In particular, 
given some specific dispersion relations, the big bounce solutions to the modified FRW equations 
can be derived. Notably, to obtain a well defined rainbow metric at the moment of the big bounce, 
we find it seems necessary to introduce a cosmological constant which depends on the energy of 
probes as well, implying that a universe with a positive cosmological constant more likely undergoes 
a big bounce at least at this phenomenological level. 
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I. INTRODUCTION 

Before a complete and fundamental quantum theory 
of gravity can be established, it has been suggested that 
a semi-classical or phenomenological theory of quantum 
gravity may play a crucial role in disclosing the nature of 
quantum gravity effects, in particular its possible impact 
on the very early universe and extremely high energy 
physics jll— flQj] - Recently such a semi-classical formalism 
named as rainbow gravity has been proposed by Magueijo 
and Smolin, which can also be viewed as an extension of 
doubly special relativity in curved spacetime [111]. In this 
formalism, one key ingredient is that there is no single 
fixed spacetime background when the quantum gravity 
effects of moving probes on geometry is taken into ac- 
count. Instead, it is replaced by a one parameter family 
of metrics which depends on the energy of these probes, 
forming a "rainbow" metric. More explicitly, suppose the 
modified dispersion relation in doubly special relativity 
has a general form 
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where two general rainbow functions / 2 (e) and g 2 {s) de- 
pend on the energy e of probes and may be expanded 
in the order of Planck length l p = v8ttG ~ l/M p . Ob- 
viously one requires that f 2 (e),g 2 (e) ~ 1 as e/M p <C 1 
. Correspondingly, it is conjectured that the usual flat 
metric should be replaced by the rainbow metric defined 
as 



ds 2 = - 
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such that the contraction between infinitesimal displace- 
ment and momenta is a linear invariant [ 1 1[ 1 1211. 



dx^pn = dte + dx° pj 



(3) 



'Electronic address: 
t Electronic address: 



y ling® ncu . cdu.cn 



The flat rainbow metric ([5J indicates that the geom- 
etry of spacetime depends on the energy of a particle 
moving on it. That is to say, even in the absence of 
gravity a moving particle with energy e would probe a 
geometry which is described by an energy-dependent set 
of orthonormal frame fields 



eo = / 1 (l P s)eo, e» = g 1 (l p e)e i , 



(4) 



where the tilde quantities are the energy independent 
frame fields probed by low energy particles. In this man- 
ner the metric can be written as 



g(s) = r 1 ab e a (s)®e b (e), 



(5) 



which can be viewed as a one family of flat metrics, char- 
acterized by £. This strategy overcomes the difficulty of 
defining the position space conjugate to the momentum 
space arising in DSR where the Lorentz symmetry is ac- 
complished by nonlinear transformations in momentum 
space (For recent progress and discussion on this issue we 
refer to (HI). 

Rainbow metric formalism can be pushed forward 
when the gravity is taken into account. First of all, a de- 
formed equivalence principle of general relativity is pro- 
posed in ll||, requiring that the free falling observers 



who make measurements with energy e will observe the 
same laws of physics as in doubly special relativity. As a 
consequence, the background spacetime is described by a 
general rainbow metric g^,i/{s). Then, through the stan- 
dard process the corresponding one-parameter family of 
connection V(s)^ which is compatible with the rainbow 
metric g^, v {s) and the curvature tensor R{s)u V \ can ' 3e 



constructed, 
equations 



leading to a family of modified Einstein's 



G^(e) = 8ttG(s)T^(s) + g^(e)A(e), 



(6) 
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where the Newton's constant G(e) and the cosmological 
constant A(e) are expected to vary with the energy e as 
well from a renormalization group point of view. 

Rainbow gravity formalism has received a lot of at- 
tention recently and other stimulating work on this for- 
malism can be found, for instance, in p^ - |25j . In par- 
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ticular, in one of the authors generalized the modi- 
fied Fricdmann- Robertson- Walker (FRW) equations orig- 
inally presented in by considering the cosmological 
evolution of probes and derived solutions in which the 
spacetime curvature has an upper bound such that the 
cosmological singularity is absent. But in [2l| the big 
bounce solution to the modified FRW equations is not 
available, which greatly depends on what kind of modi- 
fied dispersion relations we would apply. 

Another motivation of looking for non-singular bounc- 
ing solutions in rainbow gravity formalism comes from 
recent progress in cosmology. It suggests that bouncing 
universes could play a similar role as inflationary sce- 
nario in solving the well known problems in standard 
Big-Bang cosmology. Nowadays there are a lot of cos- 
mological models with a bounce solution such as the 
pre Big -Bang scenario [26jand cyclic/Ekpyrotic universe 
|27l l28f . as well as superstring cosmologyf29 |-[ 3 3|, b rane 
cosmology |34|, I35L lo op quantum cosmology (36l - l38j and 
quintom models |39l . l40j . For more details on bouncing 
universe we refer to a recent review (4l| and references 
therein. Since rainbow gravity formalism is proposed as 
a semi-classical theory in which the quantum effects of 
gravity is taken into account, we would like to ask if it 
is possible to obtain such big bounce solutions in this 
framework. Bearing this question in mind, in this pa- 
per we intend to further investigate the rainbow universe 
in a more general setting. Through some explicit model 
constructions we find the answer is affirmative. 

We organize the paper as follows. In section II, we 
derive the modified FRW equations in the framework 
of rainbow universe where the cosmological evolution of 
probes is taken into account. Then we turn to derive the 
big bounce solutions to these equations with vanishing 
cosmological constant in section III. Through specifying 
modified dispersion relations, two sorts of big bounce so- 
lutions are demonstrated. But these models contain some 
unsatisfactory features at the moment when the universe 
passing through the big bounce. In section IV we show 
that such unsatisfactory points can be overcome by in- 
troducing an effective cosmological constant. The corre- 
sponding big bounce solutions are presented as well. 



constant which is independent of spacetime coordinates. 
However, in [2~H rather than considering any specific mea- 
surement, this mechanism is generalized to study the 
semi-classical effects of particles on the background met- 
ric during a longtime process. Then the probe energy 
appearing in the rainbow metric is identified with one 
of photons or other sorts of massless particles like gravi- 
tons or inflatons which dominate the very early universe. 
Thus the evolution of energy e with the cosmological time 
need to be taken into account, denoted as e(t). As a re- 
sult the rainbow functions / and g in Eq.© depend on 
time only implicitly through the energy of particles. In 
21] we take the ansatz with g' 2 = 1. Here for our pur- 
pose / and g are chosen a priori and we will derive the 
modified FRW equations in a more general setting. 

Directly starting from the rainbow metric in Eq.([7]) 
and taking the cosmological evolution of e into account, 
we obtain the non-vanishing components of associated 
connection as 
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Next the components of Ricci tensor as well as the Ricci 
scalar can be derived as follows 
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II. MODIFIED FRW UNIVERSE IN RAINBOW 
GRAVITY FORMALISM 

The modified FRW equations in rainbow gravity for- 
malism have originally been derived in 11( . As a starting 
point, the conventional FRW metric is replaced by a rain- 
bow metric parameterized by the energy e of probes 



ds 2 = 
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dt 2 



fl 2 (e) 
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Here we only consider the spatially flat case with K = 0. 
This metric is defined in a general sense in that one is 
free to pick up arbitrary particle as a probe and for any 
specific measurement its energy s can be treated as a 



Finally, substituting all the terms above into Eq. 
obtain the modified FRW equations as 
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where H = a/ a is the Hubble parameter and p, P de- 
note the energy density and the pressure of perfect fluids 
respectively. 

Furthermore, the conservation law for the energy- 
momentum tensor reads as V(e) M T M;y = 0, where the 
covariant derivative V(e) M is energy dependent as well. 
Plugging all the components of the connection into this 
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equation leads to a modified conservation equation as 



p + 3(H-±)(p + P) = - 
9 



A(e) _ G 
8^G _ G P ' 



(14) 



Obviously the conservation equation (|14|) is not indepen- 
dent and can be derived from Eqs. (IT^|) and (fTBI . 

We now turn back to Eq. ([7]) and look at the term 
a(t)/g 2 , which depends on the energy of probes through 
the factor g~ 2 . Obviously, at low energy limit — ^ 1 , it 
is just the ordinary scale factor measured by observers. 
The cosmological singularity occurs in standard model 
when one tracks back to the origin of the universe and 
insists to take this as the physical scale factor, namely 
a(t) — > as t — > 0( or at some finite time). How- 
ever, in rainbow gravity formalism we find that as the 
energy of those particles evolving along with spacetime 
geometry increase at the very early stage of the uni- 
verse, the factor g 2 may be far from the unit and its 
effects can not be ignored. Then any possible measure- 
ment by those particles is energy dependent such that the 
actual and physical scale factor measured by a probe de- 
pends on its energy e through the function g 2 , which is 
quite different from the ordinary scale factor a(t). As 
a matter of fact, in this formalism as ait) — > 0, the 
rainbow function g~ 2 probably becomes divergent such 
that the total effective scale factor probed by particles 
may be finite, providing a mechanism of avoiding the 
singularity. Therefore we define an effective scale factor 
a e ff — a/g and subsequently the effective Hubble param- 
eter as H e ff = a e ff/a e ff = H — g/g. They are the phys- 
ical parameters probed by those matter evolving with 
space time geometry in very early universe. Through 
this paper we will focus on analyzing the solutions to the 
effective scale factor a e f / . 



III. THE BIG BOUNCE SOLUTIONS WITH 

A = 

Now we are going to look for big bounce solutions to 
the modified FRW equations (fT2)) and (fTB")) . For sim- 
plicity through this paper we will treat the Newton's 
constant G to be energy independent. Next we need 
to specify the rainbow functions /(e) and g{e). It has 
been pointed out in [2l| that the effective FRW equa- 
tions appearing in loop quantum cosmology [42l444T | can 
be heuristically derived from the framework of rainbow 
gravity once the modified dispersion relation is properly 
assigned. For instance, the effective equation for K = 
and A = in loop quantum gravity has the form 



H 2 = ^-p I 1 
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where p c ~ L 4 . They may be obtained from the side of 
rainbow gravity if we naively set that 



and e 4 ~ e 4 ~ p. This identification presumably in- 
dicates a relation between the rainbow gravity formal- 
ism and the effective theory of loop quantum gravity at 
the semi-classical level. As a result it appears that the 
big bounce solutions should be obtained in rainbow uni- 
verse straightforwardly. However, as pointed out in[2llj . 
this naive identification involves more subtleties to clarify 
from the side of doubly special relativity. The reason is 
that the modification of dispersion relations will provide 
corrections to the expectation value of statistical quan- 
tities of an ensemble as well. If a modified dispersion 
relation does not manifestly provide an upper bound on 
either the momentum or energy of a single particle, it 
would lead to a divergent density of states as the energy 
approaches to the Planck scale, which can be seen from 
the following definition (2j HBJ 
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(17) 



where prime denotes the derivative with respect to the 
energy e. To preserve the finiteness of the mean value 
of statistical quantities with the use of Eq. (IT71) , it seems 
that some certain cutoff of the energy should be intro- 
duced by hand, which obviously sounds not satisfying. 
A more detailed analysis shows that this difficulty does 
arise when one intends to derive the big bounce solution 
in rainbow universe through the ansatz (|16p . To avoid 
this weak point, we propose the following two sorts of 
modified dispersion relations and then discuss the solu- 
tions to the corresponding modified FRW equations. 



A. the case of / = g — fzjrp 

The first sort of rainbow functions we would like to 
adopt is 



f = 9 2 
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(18) 



f = (i-iyr\ 
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(16) 



Obviously in any case of f 2 = g 2 , the ordinary relation 
de = dp is preserved for massless particles such that all 
the statistical quantities will not receive corrections due 
to the modification of dispersion relations. In particular, 
the equation of state P = uip and the relation p ~ e 4 still 
hold for those particles. 

Next we identify the energy appearing in rainbow met- 
ric with the statistical mean value of all massless radia- 
tion particles, namely s ~ e. This is because we are 
concerned with the "average" effect of radiation parti- 
cles which dominate the very early universe, rather than 
picking up any specific particle from the radiation at ran- 
dom. The strategy of treating all radiation particles as an 
ensemble and considering their statistical effects has pre- 
viously been applied to investigate the thermodynamics 
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FIG. 1: The effective scale factor a e // vary with the cosmo- FIG. 2: The effective scale factor a e ff vary with the cosmo- 
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logical time when f 2 — g 2 



logical time when / = g 2 = yzji 



of black holes as well [18j, [l9j, |46j . As a result, the rainbow 
functions of / and g can be finally expressed as 



f=9 2 = 



1 



1 



1 - 



Correspondingly the modified FRW equations read as 
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(20) 
(21) 



which are exactly the effective FRW equations arising 
in loop quantum gravity. Without surprise when a; is a 
constant but not equal to —I, we have the big bounce 
solution as 



1 



P = 



and 



a eff = [{Ct) 2 



(22) 



(23) 



where C 



.3(1 



V " ^ * s eas ^ to nnc ^ ^hat tne en- 
ergy density p is bounded at the big bounce of the uni- 
verse, namely 



p —> -—j, as t 0. 



(24) 



Therefore, in rainbow universe the cosmological singular- 
ity can be avoided. For explicitness, we plot the evolu- 
tion of the effective scale factor for two special cases with 
uj = —2/3 and w = 1/3, as illustrated in Fig.l and Fig. 2 
respectively. 

From Eq. (|2"3")) . we find there is always a big bounce for 
the very early universe if u; > —1. However, some issues 
arise when one attempts to understand the behavior of 



the unverse at the moment t = 0. Though both the 
energy density and the effective scale factor are finite at 



t = 0, the rainbow functions / 
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approach 



i-ZpV 

to infinity as p — > 1/lp such that the time-like component 
of the rainbow metric vanishes, implying the degeneracy 
of the metric at that time. Such difficulty is general for 
this sort of modified dispersion relations and can not be 
avoided by inserting any adjustable parameter like f 2 = 
g 2 = -. — To overcome this unsatisfactory point, we 
consider another sort of modified dispersion relations in 
next subsection. 



B. the case of / = g = 1 + lp£ 4 

Now we attempt to take the rainbow functions / and 
9 as 

/ = g = 1 + l\e\ (25) 
Then the effective FRW equation reads as 



H eff 



BttG p 



(26) 



Combining this equation and the conservation equation 
(fT4"l) leads to the evolution equation of the energy density 
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The solution to this equation is 

(Ct) 2 + 2l p 4 - ^/(Ct) 4 + 4l p 4 (Cty 



P = 



— . 



(27) 



(28) 



where C is a constant. Again the energy density is 
bounded and approaches to jx as t — > 0. Especially, 
in this case both rainbow functions / and g at t = 
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FIG. 3: The effective scale factor a e // vary with the cosmo- 
logical time when / = g = 1 + l*p, ui — — | 



In this setting Eq.([30|) becomes 
8ttG " 



H eff ~ 



P (l - l 4 pP ) + (1 + X P )(1 - l 4 pP ) 



A 



8nG 



Combining Eqs. (fT4]) and (|3T1) . we have 



(31) 



= -M, 



1 + 7^ 0. Introducing new variables 



p^-n + x S)p 2 + [i + - mp + 

(32) 

where M = 3(1 + u)^^ S J£+ XA ■ Here, we require 

8ttG 
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PA 



A 
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Eq. (|32|) can be changed into 



are finite such that the rainbow metric is well defined 
even at the moment of passing through the big bounce. 
Furthermore, we may obtain the effective scale factor as 



{Ctf 



2/ p 



V(ct) 4 + ubiety 



21 1, 



(29) 
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It is also bounded and approaches to as t — » 0, 

thus a big bounce occurs whenever For instance, 

we illustrate its evolution with u> = —2/3 in Fig. 3. From 
this figure, however, we see that the evolution at the 
moment of big bounce is not smooth. This comes from 
the fact that the solution (|29|) contains a square root such 
that a e ff is not continuous at t = 0. 

In next section, we intend to demonstrate that all the 
unsatisfactory points appearing in this section can be 
overcome by introducing a non-zero cosmological con- 
stant. 



IV. THE BIG BOUNCE IN THE PRESENCE OF 
THE COSMOLOGICAL CONSTANT 

In the framework of rainbow gravity, it is reason- 
able to expect that both the Newton's constant and the 
cosmological constant are dependent on the energy of 
probes [llj. Here we only consider the possible modifica- 
tion of the cosmological constant A(e). Introduce a third 
rainbow function h such that A(e) = h 2 (e)A where A is a 
bare cosmological constant, the modified FRW equation 
(fT2"|) can be written as 



H eff 



8nGp h 2 (e)A 
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(30) 



To obtain the big bounce solutions and compare them 
with the results in previous section, we fix the functions 



/ and g as / 
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and assume h 2 = 1 + \p. 



^p A x 2 + [1 + (A - lj) PA \x - /4(l + Ap A ) 



= M. (34) 



Next we present solutions to this equation in accord with 
the sign of the cosmological constant A. Here for conve- 
nience, we define T as 

T = PA x 2 + [1 + (A - l A p ) PA ]x - + Xp A ). (35) 



A. The case of A > 

When A > 0, we obtain the solution to Eq. (|34|) as 

MJp^t = ln[2p A .T + [1 + (A - 11) P a] + 2VPA 

/T]-lniV, (36) 



where to make the big bounce occur at t = we have 
introduced an integral constant N > which is defined 



as 



[1 + (A - l 4 p ) PA ] 2 + 4Z> A (1 + A PA ) > 0. (37) 

If we consider A > — Eq. (|37|) has been satisfied. It 
is very useful to rewrite the solution in Eq. (|3"o]) ~ 



as 



= - — exp (MJ~p~At + In N) 
4pA 

+ — N 2 exp[-(M^/pXt + lnJV)l 
4pA 

~ [1 + (A ~^ )/9A] 2^' 



(38) 



It is easy to show that x has the minimal value at t = 0, 



Xmin = ^-[N-l-{X-l 4 )p A ] 

2pA p 

IP + Apa) 
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FIG. 4: The effective scale factor a e // varies with the cosmo- 
logical time when A > 0. 



FIG. 5: The effective scale factor a e ff varies with the cosmo- 
logical time when A < 0. 



If the parameter A satisfies the condition 



A > 



1 

PA 



(40) 



then we find the energy density p < 4- at the big bounce. 

As a matter of fact, in rainbow gravity formalism it is rea- 
sonable to expect the parameter to be an order of |A| ~ Z 4 , 
the condition (|40|) can be easily satisfied. Therefore, in 
the presence of a cosmological constant the rainbow func- 
tion / will never diverge such that the rainbow metric is 
always well defined even at the big bounce. Moreover, 
the effective scale factor evolves as 

a e ff = j^exp(MVM« + lniV) 

+ — N 2 exp [- {MJpZt + In N)} 

4/9A 



[1 + (A-Z> A ] 



l + Ap A 
1 "| 3(l + w) 



(41) 



We plot its evolution in Fig. 4, from which we find the big 
bounce occurs at t = 0. 



B. The case of A < 

When pa < 0, we obtain the solution by integrating 
Eq.®, 

ft . /-2p A z-[i+(A-j> A n 

M t = \ arcsm < > 

V PA I N | 



7T f-i 

2 V pi' 



(42) 
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2 V PA 

that the big bounce shifts to t = 0. Similarly, we rewrite 



where 5 A / — is introduced as an integral constant such 
2 y pa ° 



this solution as 
1 



x = -^- K {Nsm(MV=£Xt - E) + [l + (A - ll) PK ]] . 

(43) 

Obviously x has the minimal value at t = 0, 



Zrmn = — [JV - 1 - (A - l*) PA ] 

2pA ^ 
g(l + A PA ) 



l + (A-I*)p A " 
If the parameter A is constrained as 



— <A<-— + /* 

PA PA P 



(44) 



(45) 



then it is guaranteed that the energy density will never 
exceed the Planck energy density, namely p < l~ 4 and 
the rainbow metric is always well defined. 
The solution of the effective scale factor is 



1 



{- — [Nsin(My/=pJt-- 



l + Ap A 



+ (l + (A-/> A )]}^. 



(46) 



We illustrate its evolution with the cosmological time in 
Fig. 5. This is an oscillation solution containing many big 
bounces. 

Similar big bounce solutions can be found in the 
context of brane scenarios and effective loop quantum 
cosmology |37ll38| . where they considered a p 2 -type mod- 
ification to the Friedmann equation due to the brane ef- 
fects or discreteness of quantum geometry. Here we ob- 
tain the big bounce solutions by considering the quantum 
gravity effects of probes. By appropriately choosing the 
rainbow functions / 2 , g 2 and h 2 we find that p 2 -type 
modification to the Friedmann equation can be achieved 
in the framework of rainbow universe as well. But in con- 
trast to brane cosmology or loop quantum gravity where 
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the modification emerges from the brane effects or quan- 
tum geometry such that the coefficients in all correction 
terms are uniquely fixed, here we treat the rainbow func- 
tions f 2 , g 2 and h 2 as independent modification quanti- 
ties. It is also interesting to notice that in loop quantum 
gravity the lattice regularization need to be refined in 
order to assure the big bounce occurring at the Planck 
scale |42h44| . While in rainbow gravity formalism, we find 
that the magnitude of the effective scale factor at the big 
bounce is adjusted by the parameter A in rainbow func- 
tion h which is assumed to be the order of l*. We may 
also compare our solutions with those appearing in low 
energy effective theory of string theory, for instance [33[ , 
where the existence of no-singular solutions is quite rele- 
vant to the dilatonic coupling constant. In our cases the 
specific form of rainbow functions plays a similar role, 
but the presence of the cosmological constant makes it 
possible to obtain bouncing solutions as well as oscilla- 
tory solutions. 

V. CONCLUSIONS 

In this paper we have constructed a phenomenological 
model for the very early universe in the framework of 
rainbow gravity, in which the quantum gravity effects of 
probes on spacetime background is taken into account at 
the semi-classical level. Starting from the rainbow Ein- 
stein's equation, we derived the modified FRW equa- 
tions by considering the cosmological evolution of parti- 
cles. Furthermore, given appropriate rainbow functions 
we find that the big bounce solutions to the modified 
FRW equations can be obtained. Notably, to obtain a 
well defined rainbow metric at the moment of the big 
bounce, we find it seems necessary to introduce a cosmo- 
logical constant which depends on the energy of probes 
as well. This is a quite interesting result. It implies that 
a universe with a positive cosmological constant more 
likely undergoes a big bounce at least at this phenomeno- 
logical level. A complete understanding on the behavior 
of the universe at the big bounce probably calls for a 
complete quantum theory of gravity. However, our semi- 
classical analysis presented here has shed light on this is- 
sue and indicates that such quantum gravity effects can 
contribute important corrections to the classical general 
relativity so as to provide a more reasonable picture on 
the origin of the universe. The validity of semi-classical 
analysis has also been strengthened by recent progress 
in loop quantum gravity in which the exact solutions 
to quantum equations are remarkably well approximated 
by a naive extrapolation of some semi-classical effective 
equations [42M-44I ] . 



Recently bouncing cosmology has received a lot of at- 
tentions as the astronomical observation has released the 
cosmological data with much higher precisions, which 
provide strong evidences for a nearly scale- invariant spec- 
trum of primordial density fluctuations [47|. First of all, 
theoretical investigation indicates that it is possible to 
generate such a spectrum during the contracting phase in 
various bouncing models, for instance see references [48l— 
[Hot ]. Secondly, in contrast to the slow-roll inflationary 
scenario in which the non-Gaussian contribution to the 
density fluctuation spectrum is greatly suppressed, some 
bouncing models require that the universe experience a 
contracting phase before the hot Big Bang expansion 
such that a large non-Gaussianity can be predicted [5ll— 
The forthcoming observations maybe distinguish 
and rule out various models of the origin of the universe. 
The non-singular bounce solutions that we have obtained 
in this paper make it possible to investigate the perturba- 
tion theory of the bouncing cosmology in rainbow gravity 
formalism. Our investigation along this direction is un- 
der progress. 

Through the paper we only investigate some special 
cases in rainbow gravity with specified functions of /(e), 
g[e) and h(e). However, we would like to point out that 
a large class of rainbow functions may lead to bouncing 
solutions if they could provide a /9 m -type modification to 
the Friedmann equation where m is not necessarily fixed 
to be square as in references. Such contributions would 
suppress the effective Hubble parameter until it vanishes 
at the big bounce. Of course, we expect the further tests 
relevant to Lorentz symmetry would tell us which is the 
proper one among all the possible modified dispersion 
relations. 

It is completely possible to extend the scheme pre- 
sented here to more general cases, for instance when the 
spatial metric is non-flat and the Newton's constant is 
also energy dependent. 
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